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We present two approaches to the dynamics of a quench- induced phase transition in the quantum 
Ising model. One follows the standard treatement of thermodynamic second order phase transitions 
but applies it to the quantum phase transitions. The other approach is quantum, and uses Landau- 
Zener formula for transition probabilities in avoided level crossings. We show that predictions of 
the two approaches of how the density of defects scales with the quench rate are compatible, and 
discuss the ensuing insights into the dynamics of quantum phase transitions. 
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Studies of phase transitions traditionally focussed on 
equilibrium scalings of various properties near the crit- 
ical point. The first major exception was an attem 
to model the physics of the early Universe: Kibble 
noted that cosmological phase transitions in a variety of 
field theoretic models lead to formation of topological 
defects (such as monopoles or cosmic strings) which may 
have observable consequences. One of us then pointed 
out that analogues of cosmological phase transitions 
can be studied in the laboratory. In such experiments 
the equilibrium critical scalings predict various aspects 
of the non- equilibrium dynamics of symmetry breaking, 
including the density of residual topological defects 

These ideas led to the Kibble- Zurek mechanism (KZM), 
a theory of defect formation that uses the critical scalings 
of the relaxation time and of the healing length to deduce 
size (£) of domains that choose the same "broken sym- 
metry vacuum" 13-01 When the broken symmetry phase 
permits their existence, KZM predicts defects will appear 
with density of about one defect unit (e.g., one monopolc 
or a £-sized section of a string) per £-sizcd domain. This 
KZM prediction has been tested, extended and refined 
with the help of numerical simulations 0, |(| , and verified 
in a variety of increasingly sophisticated and reliable ex- 
periments in liquid crystals [J El, superfluids lit [Tel. Ill| . 
superconductors [13, LUi > and other systems [lj| . 

A majority of the experimental data agree with KZM. 
One notable exception is the case of supcrfhiid 4Hc, 
where initial reports of KZM vortices being detected 
were retracted [lfj after it turned out that stirring had 
inadvertently induced vorticity. In view of various un- 
certainties, it is still not clear whether 4He experiments 
are at odds with the numerics-assisted KZM predictions. 
Regardless, KZM provides a theory of the dynamics of 
second order phase transitions ranging from low tem- 
perature Bose-Einstein condensation to grand unification 
scales encountered in particle physics and cosmology. 

In this paper we consider a barely explored problem: 
the dynamics of quantum phase transitions. Quantum 



many-body systems (e.g., Bose gases) can undergo ther- 
modynamic phase transformation (such as Bose-Einstein 
condensation that follows evaporative cooling). KZM 
theory, developed to deal with thermodynamic phase 
transitions, applies in this case directly, even though the 
dynamics of Bose condensation is explicitly quantum 1 161 . 

On the other hand, a quantum phase transition |l7j 
(e.g., the Mott insulator-supcrfmid transition of bosons 
in a periodic lattice) is a change in the character of a 
system's ground state which occurs as some parameter of 
its Hamiltonian passes critical value. For instance, low- 
ering of the amplitude of the optical lattice induces Mott 
transition. Unlike thermodynamic transitions, quantum 
phase transitions involve only reversible unitary dynam- 
ics. Therefore, scaling arguments that work in thermody- 
namic transitions (where the order parameter is damped) 
may not be valid in the quantum case (but see [l8lll9| ). 
Yet, dynamics of quantum phase transitions is interesting 
in its own right and has applications in quantum infor- 
mation processing 0, [2(| ■ 

We will study a quench-induced transition in the quan- 
tum Ising model. This model is regarded by Sachdev 01 
as one of two canonical quantum phase transitions. It de- 
scribes a chain of spins with the Hamiltonian: 
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Here erf, erf are Pauli operators, W is the Ising coupling, 
and J(t) is proportional to the strength of an external 
field that attempts to align spins with the x-axis. 

The phase transition from the high-field state (all spins 
aligned with x, i.e., | —►,—►,...,—►)) to the low-field 
ground state manifold - spanned by | TjT> ■■•>!) and 
I I, I, . . . I), and doubly degenerate in the large N limit 
- takes place when J(t) = W. Thus, relative coupling: 



<(t) = J(t)/W - 1 



(2) 



is expected to play a role of the relative temperature 
— 1) in system's behavior near the critical point Tc- 
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Indeed, all relevant properties depend on the size of the 
gap (A) between the ground state and the first excited 
state. As N — ► oo, the gap is: 

A = 2\J(t) -W\ = 2W\e(t)\ . (3) 

and sets an energy scale reflected in the relaxation time 

t = h/A = h/2W\e(t)\ = T /\e(t)\ . (4) 

Divergence of r near the critical point is the critical slow- 
ing down familiar from thermodynamic phase transitions. 

The healing length is given by the product of the speed 
of sound (c) and the relaxation time: 



Z = 2Wa/A(t) = a/\e(t)\=Z /\e(t)\ 



(5) 



where c = 2Wa/h (see [lTj), and a is the distance be- 
tween spins. The divergence of £ near the critical point 
is analogous to critical opalescence. 

The scaling of r and £ suggests estimating the size of 
broken symmetry domains (i.e., regions of aligned spins) 
using the same approach that worked in thermodynamic 
transitions 0, Q : Near the critical point "reflexes" of the 
system (measured by the relaxation time r) deteriorate, 
until - at the critical point, where r = oo - system cannot 
react at all. Yet, early in the quench r is still small, 
and its state is still able to adjust to variations of the 
external parameter (e.g., J or T). This suggests splitting 
the quench into the near-critical impulse regime and the 
quasi-adiabatic regime far from the critical point. Such 
split is the essence of KZM. 

The instant (t) when behavior of the system changes 
from adiabatic to impulse is of key importance. This hap- 
pens when its reaction time (given by Eq. is the same 
as the timescale on which its Hamiltonian is changed. To 
calculate t, we assume that the external bias field changes 
linearly with time, so that e(i) = £/tq. As the relative 
coupling changes on a timescale e(t)/e(t) = t, the switch 
between adiabatic and impulse regimes occurs at the in- 
stants ±i when relaxation time is equal to t, 

r(£)= r /|e(£)| = r T Q /t = e(£)/e(f)=f (6) 

which yields 

i = = \jrqhj2W . (7) 

Typically, these two instants (±£) separate evolution 
into three regimes. Initially, for t < —t, the system's 
state will adjust to the decreasing J(t). However, at 
t = —t (before the critical point) this tracking of the 
instantaneous ground state of H will cease. Evolution 
will re-start only at +t (after the critical point), with an 
initial state similar to the one "frozen out" at —t. 

In thermodynamic phase transitions fluctuations of the 
order parameter at t give rise to domains of size £ given 
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FIG. 1: Density of kinks (y) in quantum Ising model after a 
quench that starts in a ground state at J = 5W and ends at 
J — 0, plotted versus the dimensionless quench rate to/tq = 
hv/AW 2 for N = 50,60,70,80,90,100 (solid lines; bottom 
to top). The simulation data are consistent with the scaling 
vkzm ~ \J to I tq , Eq. 1)10(1 . predicted by the Kibble-Zurek 
mechanism. (See [l9|l for details on the numerical method). 
Agreement improves with N; for 100 spins, a fit gives v ~ 
Tq (dashed line). As in the classical case @ Eq. JTSJ is 
an overestimate; the best fit is v ~ 0.16 vkzm- 



by the healing length at —t. Using the relative coupling 
e at t we similarly calculate for the quantum Ising model: 



e{t) =t/T Q = J to/tq 



(8) 
(9) 



£ = Co/e = £o\ tq/tq = aJ2WT Q /h 



Note that this scaling differs from the £ = £o/Ve 
£o( r Q/ T o) 1 predicted by non-relativistic mean-field the- 
ories for second order phase transitions 0, El ■ 

Following KZM, we now expect appearance of 0(1) 
defects per £. Their density should be therefore: 



vkzm 



a/i = ^Jh/2W 



T Q- 



(10) 



This is only an estimate. Simulations of classical second 
order transitions yield defect densities that scale with tq 
as predicted by KZM, but that can be lower by about an 
order of magnitude than a/£: Defects can be separated 
not by £ but by approximately 10 | (see @). 

This KZM paradigm should not be uncritically applied 
to quantum phase transitions. Above all, thermodynamic 
fluctuations are 'real'. If they survive from before -t, they 
can tip the balance at +t, breaking symmetry right after 
the transition. It is hard to make an analogous argument 
for a quantum case equally convincing. Quantum fluctu- 
ations, exist, but they are virtual, so it is not obvious 
that they will have a similar symmetry breaking effect 
on the post-transition state. On a more prosaic note [as 
we shall see in Fig. Gla)] the relevant gap (i.e., gap be- 
tween the ground state and the first accessible excited 
state) is not the symmetric A, Eq. ©. Rather, its slope 
is twice as large on the approaching side. 

Nevertheless, Fig. ^ shows that number of kinks per 
spin - the residual kink density created by quenching the 
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FIG. 2: (a) Energies of lowest excitations (see [Tfll p for N = 20. The energies of the ground state and the first accessible excited 
state are the lowest (horizontal) solid line and the second lowest solid line, (b) (i) Quench time tq/tq = 4W 2 /hv that yields 
/ of 99%, and (ii) the fidelity for a fixed tq = 200h/W = 400to as a function of the number of spins N in the quantum Ising 
chain. A power-law fit to the data corresponding to tq 99% gives a power of 1.93 (LZF yields 2, as would KZM). The best fit for 
the fidelity with Landau-Zener dependence / = 1 — exp {— oWtq/HN 2 } yields a ~ 59, compared to theoretical a — 2-k' a ~ 62 
based on the steeper slope in the marked energy level above (a = 4ty 3 ~ 124 when the more shallow slope in Fig. HJa) consistent 
with Eq. (3) is taken), (c) Upper (dashed lines) and lower (solid lines) bounds on fidelity as a function of to/tq = hv/4W 2 for 
N = 90 (i), N = 70 (ii), N = 50 (iii), N = 30 (iv). Fits to LZF (dotted lines) lie between these bounds, and for / > 0.6 give 
a~ 59 (i), 59 (ii), 57 (iii), 54 (iv). 



quantum Ising model scales approximately as ~ 
just as predicted by Eq. iJTUJ. This holds throughout the 
region of KZM's validity, i.e. where e is much less than 1 
(so that the quench is quasi - adiabatic early on and at the 
end, but 'impulse - like' near the critical point, and thus 
at least one defect is expected). The prefactor ~ 0.16 
[0.12 if the steeper slope on the approach in Fig.EI a ) is 
taken] is also not far from previous experience |5j . 

Kibblc-Zurek mechanism works in a quantum transi- 
tion! Yet, in view of doubts about quantum fluctuations, 
an explicitly quantum treatment would be reassuring. 

As N — > oo, the gap (A) - a salient feature of the 
quantum Ising model - disappears at the critical point 
in accord with Eq. |J5J. When N < oo, this critical gap is 
small, but finite [see Fig. 2(a)]. This is of key importance 
for the remainder of our paper. Instead of calculating the 
density of defects in an infinite system, we shall compute 
size (JV) of the largest spin chain likely to remain defect 
free (i.e., in a ground state) after a quench, as a function 
of quench timescale tq. For defect probability of ~50%, 
inverse of N is an estimate of defect density. 

Excited eigenstates of H, Eq. (JJJ, on the broken sym- 
metry (W > J) side of the transition describe states of 
the spin chain in which the direction of symmetry break- 
ing varies along the chain once, twice, etc. [T^ |. Thus, 
they represent states containing one, two, etc. "kinks" . 
The behavior of the energies of lowest excitations of H 
in the vicinity of the critical point, [Fig. EI a,)], suggests 
avoided level crossing. Hence, it appears that phase 
transition dynamics in the quantum Ising model can be 
treated using the Landau-Zener formula pl|. or LZF. 
LZF gives the probability of exciting a system driven 
through an avoided level crossing: 



Here, A is the minimum energy gap between the two 
levels, and v is the quench velocity. That is, far away 
from the "point of the nearest approach", v = A. 

Using LZF we can compute the average size N of a 
spin chain that is likely to remain in the ground state 
throughout the quench. In the adiabatic limit (v ~ 0), 
Eq. (|11|) predicts that the system will stay in the same 
energy eigenstate (i.e., the probability of switching levels 
will be vanishingly small). To quantify this Ref. [19j uses 
the fidelity, / = \{iI)actual\^ground)\ 2 , which gives 
the probability that no defects will be produced. From 
LZF it follows that pchange - cx p(-§^|) - 1 ~ / • 
Thus, the rate of a nearly defect - free quench (resulting 
in defects with probability 1-/ <C 1) is bounded: 



ttA 2 



2h\ ln(l-/)| 



(12) 



Below wc will express v using quench time tq and W; 
v = |A| = 2J(t) = 2W/t q . 

The lowest excited states are inaccessible - they have 
a different parity than the ground state, and H conserves 
parity. The first accessible level has one kink for J = 0. It 
gets to within A = AnW/N for N 3> 1 above the ground 
state. With these ingredients using Eq. (|12|) wc obtain: 



IAI = 



2W 7r(47rW7A0 5 



tq ~ 2S|ln(l-/)| 



(13) 



It relates the size N of a defect-free chain to quench time: 



N < 2tt, 



itWtq 



fi|ln(l-/)| 



27rl^, 



2tt 



HMW)I 



(14) 



Pchange ^ cxp( 



ttA 2 
'2h\v\ 



(11) 



Figs. 2(b) and 2(c) show that LZF provides a good fit for 
/ > 0.5. This in not completely unexpected: Damski j2^] 
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FIG. 3: Number of kinks in chains (iV=50,60, 70,80,90, 100, 
spins, bottom to top) after a quench versus the quench rate 
r a/ r Q = hv/iW 2 . Both the scaling v ~ I/^/tq predicted by 
KZM (red lines), Eq. JTUJ, and the LZF estimate 1 - / (blue 
lines) are valid where expected. The red lines are linear fits in 
the range (0.025,0.25) yielding slopes between 0.66 and 0.58. 
The blue lines are the fit results from Fig. |2|c). Numerical 
data (obtained using the method developed in Hi! ') include 
these used in Fig. 1. However, we now go beyond the expected 
range of validity of KZM: For quenches slow enough to create 
'less than a kink' LZF provides reliable predictions, while very 
fast quenches are "all impulse", and - according to KZM - 
the number of kinks should saturate, as is indeed seen. 



recently proposed a 'KZM approximation to LZF' in an 
insightful paper. However, even upon closer inspection 
(work in progress |23j L such agreement between LZF in 
the setting which is not a standard avoided level crossing 
and the numerics might be still somewhat surprising. 

When we compare the KZM and LZF predictions for 
defect density, we find: 



1 1 

iV 2tt 



2|m(l-/)| 



x vkzm 



(15) 



The two estimates of defect density exhibit the same scal- 
ing with the quench rate and with the parameters of H, 
Eq. JTJ. However, LZF predicts fewer defects than the 
"raw KZM estimate" {vlzf — 0.14 x vkzm when / is set 
- somewhat arbitrarily - to 0.5). This is no surprise; nu- 
merical simulations, experiments and analytic solutions 
to specific models, have shown that Eqs. © and l|l(J|l 
provide correct scalings, but tend to overestimate densi- 
ties (see e.g. 0, Q). Fig. 1 confirms that this is also 
true also for the quantum Ising model. 

We note that while vkzm and vlzf are closely related, 
they answer somewhat different questions. In particular, 
vkzm does not depend on /. However, when less than 
one defect is expected in a chain, the number of defects 
is ~ 1 — / and can be computed using LZF. Fig. 3 shows 
that LZF and KZM complement each other in this case, 
and jointly cover a wide range of quench rates. 

We found that a quantum analogue of KZM, based on 
critical scalings, predicts the results of numerical simula- 
tions. As expected, KZM scaling holds when e < 1 - i.e., 



when the quench starts and ends in the adiabatic regime, 
but becomes impulse - like near the critical point. When 
the quench is so slow that it never acquires impulse - like 
character, LZF is accurate. We conclude that the two ap- 
proaches work well in complementary regimes of quench 
rates, and predict the same scaling for the size of broken 
symmetry domains with quench time. 
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